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ABSTRACT 



The standard formula for the rotation measure, RM, which determines the 
position angle, if) = RMA 2 , due to Faraday rotation, includes contributions only 
from the portions of the ray path where the natural modes of the plasma are 
circularly polarized. In small regions of the ray path where the projection of the 
magnetic field on the ray path reverses sign (called QT regions) the modes are 
nearly linearly polarized. The neglect of QT regions in estimating RM is not 
well justified at frequencies below a transition frequency where mode coupling 
changes from strong to weak. By integrating the polarization transfer equation 
across a QT region in the latter limit, I estimate the additional contribution Aif> 



needed to correct this omission. In contrast with a result proposed by Broderick 



& Blandford (2010), is small and probably unobservable. I identify a new 



source of circular polarization, due to mode coupling in an asymmetric QT region. 
I also identify a new circular-polarization-dependent correction to the dispersion 
measure at low frequencies. 



Subject headings: magnetic field - circular polarization - radiative transfer 



-3- 
1. Introduction 

Faraday rotation is the rotation of the position angle, ip, of linear polarization 
of radiation passing through a magnetized plasma whose natural modes are circularly 
polarized. The handedness of the polarization of the natural modes is an odd function of 
cos# = ay | cos#|, where 9 is the angle between the ray path and the magnetic field. The 
sense of Faraday rotation reverses when the sign ay reverses. The net amount of Faraday 
rotation is proportional to the square of the wavelength, A, and is written RMA 2 , defining 
the rotation measure, RM. The conventional formula for RM, cf. ([I]) below, has positive and 
negative contributions from the portions of the ray path with ay = ±1, respectively. This 
standard formula neglects the contributions of the small regions around the points where 
ay changes sign, and 9 passes through it/2. As such a magnetic- field reversal is approached 
the polarization of the natural modes becomes increasingly elliptical, as | cos#| decreases, 
becoming linear at 9 = n/2. Mode coupling needs to be taken into account in such regions, 
and there are two limiting cases. In the limit of strong mode coupling, the width of the 
region where the modes are not circularly polarized is negligible, and it can be treated as 
a sharp boundary. In the opposite limit of weak mode coupling the components in the 
two modes remaining in those modes, evolving from circular through elliptical to linear at 
9 = 7r/2, through elliptical back to circular in the opposite sense. The standard formula for 
RM effectively assumes strong mode coupling. The strength of mode coupling decreases 
with decreasing frequency, and below a transition frequency, Ut, mode coupling becomes 
weak. The motivation for the present investigation arose in connection with potentially 
observable effects on RM for v < v^. A new result that emerges from the investigation is 
that partial mode coupling can result in incident linearly polarized radiation emerging with 
a circularly polarized component. 

There is an extensive literature on mode coupling in connection with solar radiophysics, 
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e.g., as summarized is several monographs (Budden 1966 Zheleznyakov 1970 Melrose 



1980 Benz 2002). Although the jargon used is outdated and confusing, I use aspects of 
it here. Specifically, the approximations in which the natural modes are nearly circularly 
polarized and nearly linearly polarized are referred to as the QL approximation and the 
QT approximation, respectively. The region about the point on a magnetic field reversal 
where the modes satisfy the QT approximation is referred to as a QT region; the center 
of the QT region is at 9 = n/2. In the solar case the radiation of interest is initially 
circularly polarized. In the strong coupling limit an initially circularly polarized component 
preserves its polarization on crossing a QT region, corresponding to flipping from one 
mode to the other. In the weak coupling limit, an incident circularly polarized component 
follows the polarization of the mode it is in, emerging from the QT region with the opposite 
handedness. Intermediate mode coupling causes partial conversion of circular into linear 



polarization, leading to a net depolarization (Zheleznyakov & Zlotnik 1964). 



In the case of Faraday rotation, the initial polarization is linear and corresponds 



to an equal mixture of two circularly polarized modes. In a recent paper, Broderick fc 



Blandford ( 2010[ ) (hereafter BB) discussed the contribution of QT regions to Faraday 



rotation in the limit of weak mode coupling, which BB referred to as the "super-adiabatic" 
(SA) regime. BB argued that a different formula for RM applies in the SA limit, cf. ^ 
below, and they suggested that observation of this effect would provide new information 
on the properties of the interstellar medium (ISM). Part of BB's justification for their new 
formula (|2j) is that a "reflection" in the sense if) — > —if> occurs at each QT region in the 
SA limit. This is not correct. BB's assumption that the position angle of the radiation 
incident on the QT region can be specified arbitrarily, by if) = ipi n say, turns out to be 
incompatible with the constraint that the incident radiation be an equal mixture of the 
two modes. This constraint determines if> at the center of the QT region, and the transfer 
equation determines if) everywhere else along the ray path (in the SA limit). One is not 
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free to specify ip- m arbitrarily at any point identified as the boundary of the QT region: it is 
predetermined by the value at the center of the QT region and the details of the model. A 
more qualitative argument suggests that ip — > —ip cannot be correct. Suppose there is a net 
total of N complete rotations along the ray path when the first QT region is encountered, 
implying ip = ip in + ttN, where for the sake of discussion I assume that ip- m < it can be 
specified. The QT region has no "knowledge" of the number N, and cannot reverse this 
through ip- m + ttN — » — ip- m — ttN, as BB's ip — > —tp requires. The actual evolution of ip 
through a QT region (in the SA limit) is described in detail below. 

My purpose in this paper is to give a critical review of the effect of a QT region on 
Faraday rotation, and its generalization to the case where the modes are not circularly 
polarized. A specific motivation is to justify the foregoing comments on BB's treatment. 
More generally, I consider the possibility that a QT region can induce a net circular 
polarization, and find that this is indeed possible. Most of my detailed discussion is 
restricted to the limit of weak mode coupling (SA limit), which is seemingly simple: the 
radiation initially in one mode remains in that mode and the polarization of that mode at 
any point is determined by the local parameters. An important point that distinguishes 
the case of Faraday rotation from other examples of mode coupling is the constraint that 
the initial radiation is an equal mixture of the two (oppositely circularly polarized) modes. 
This constraint leads to the surprising (to me) implication that tp- m cannot be arbitrarily 
specified. I describe the transfer of polarized radiation using a matrix equation for the 



Stokes parameters, I,Q,U,V (Melrose & McPhedran 1991). This formalism does not 



involve separating into the two natural modes, and can be used without any reference 
to modes or mode coupling. The properties of the modes are implicit and can be made 
explicit by identifying the four eigenvalues of the transfer matrix, and constructing the four 
eigenfunctions (combinations of /, Q, U, V). Two of the eigenvalues are zero, and the other 
two correspond to the natural wave modes. A zero eigenvalue implies a conserved quantity, 
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and / = constant is one of these. The other conserved quantity is determined, for Faraday 
rotation across a QT region, by the constraint that the radiation be an equal mixture of 
the two modes. To treat Faraday rotation in the SA limit I derive a differential equation 
for tj) = | arctan({7/(3) subject to this constraint. 

Some background results are summarized in Section [2] polarization transfer is discussed 
in Section [3j the existence of circular polarization due to mode coupling at a QT region is 
pointed out in Section |4j the possibility of observing the effects identified is discussed in 
Section |5j and the conclusions are summarized in Section [6] 

2. Phases, modes and the Faraday angle 

In this section I write down results relating to Faraday and mode coupling that are 
needed in the subsequent discussion. 



2.1. Standard formula for RM 

The standard formula for RM can be written the form, BB's Eqn (15), 

RM = C J n e B ■ d£, (1) 

where B • di is proportional to cos 6, and where C depends on fundamental constants. The 
integral in ([!]) is along the ray path through the QL regions. The contribution of the QT 
regions is neglected. 

Although the derivation of ([I]) is well known, there are several points worth emphasizing. 
Formula Q may be derived by separating initially linearly polarized radiation into two 
equal circularly polarized components, corresponding to the two natural wave modes, 
labeled o and x say. The difference in refractive index, n Q — n x > 0, causes the two 
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components to get out of phase at a rate d(<p a — (fi^/di = Ak = (n Q — n x )u/c per unit 
distance, £, along the ray path. An important point is that the handedness of the o and x 
modes depends on the sign oy = cos#/| cos#|, and that Faraday rotation results from the 
difference in the refractive indices of the right (R) and left (L) handed natural modes, 
which is tir — ul = oy(n — n x ). Within a QL region, the position angle rotates at the 
rate dip/d£ = |oyAA;, which integrates to give ([!]). The sign ay is opposite in successive 
QL regions, with the reversal occurring at the center of the QT region that separates the 
successive QL regions. 



2.2. BB's formula for RM in the SA limit 

The new formula proposed by BB to take account of the effect of the QT regions in the 
SA limit is, BB's Eqn (14), 

RM SA = C J n e \B-d£\. (2) 

As in Q the integral is along the ray path through the QL regions. The contributions of 
the QT regions do not appear explicitly in (J2|. 

I give the following physical interpretation of (|2j). If the phase difference between the 
o and x modes is integrated along the ray path through the QL regions, neglecting the QT 
regions, the accumulated phase difference is RM S aA 2 /2, with RM S a given by (J2|. Suppose 
that the SA limit applies at every QT region along the ray path. Then the two components 
defined by the initial separation into modes remain in these modes everywhere along the 
ray path, including every QT region. In this case the accumulated phase difference between 
the o mode and x mode components corresponds to the actual accumulated phase difference 
that develops between the two components defined by the initial separation at the source. 
It is tempting to think that in the SA limit this accumulated phase difference is relevant 
to Faraday rotation, as ^ appears to imply. However, Faraday rotation is determined by 
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the accumulated phase difference between the R and L modes, and not by that between 
the o and x modes. What happens in the QT regions does not affect the fact that the 
contribution of the QL regions is given by Q, provided only that the contribution of the 
QT regions is ignored. BB's argument for ^ in the SA is that the contribution of each QT 
region is to cause ip — > —ip, which effectively converts (IT]) into fl2j). 



2.3. Mode coupling 



Traditionally, the strength of mode coupling is described by a coupling coefficient, Q x , 
such that mode coupling is strong for Q x 3> 1 and weak for Q x <C 1. The theory of mode 



coupling leads to (Cohen 1960) 



Q, 



V 



7T V 2 p V%L 



1/4 



(3) 



where L is a characteristic length over which 9 changes, v v is the plasma frequency and 
Vb is the electron cyclotron frequency. The strong dependence on v/vt suggests that the 
transition from strong to weak mode coupling occurs rapidly as a function of decreasing v. 



3. Transfer of polarized radiation 

In this section I present the matrix method for treating the transfer of polarized 
radiation in the SA limit and apply it to Faraday rotation. 



3.1. Transfer of the Stokes vectors 

The transfer of polarized radiation between the source and the observer can be 
described by the transfer equation for the Stokes parameters, /, Q, U, V. In the present 
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case, J is a constant. The transfer equation for the polarization can be written in the 



matrix form (Melrose & McPhedran 1991) 



u 

\V J 



( -pv 

Pv 

\ ~PU PQ 



PU 
-PQ 

/ 



U 

\V J 



(4) 



The parameters pQ,pu,py are derived from the (cold plasma) dielectric tensor for the 
plasma. The transfer equation Q is derived for constant Pq,Pu,Pv, an d it remains valid 
provided that any variation in pQ,pu,pv is slow compared to the rate, Ak, at which the 
modes get out of phase. 

An alternative way of writing Q involves interpreting S = (Q, U, V) and 
P = (PQ, Pu, Pv) as vectors in a 3-dimensional vector space. Then Q becomes 
dS/d£ = p x S. This alternative form implies two conserved quantities, S 2 = Q 2 + U 2 + V 2 
and p ■ S = pqQ + pjjU + py V. In a slowly varying anisotropic medium, p varies slowly in 
both magnitude and direction. The adiabatic invariants are then S 2 // 2 and p- S/pI, with 
p = \p\ = Ak. (I varies as the refractive index varies, but this is insignificant here.) The 
invariant p ■ S/p has important implications for Faraday rotation in the SA limit. 

No separation into modes is made in deriving (|4]), and no separation into modes is 
needed to use it to describe Faraday rotation. In principle, given a model for the medium, 
one can integrate ^ along the ray path from the source to the observer. Some examples 
of numerical solutions of ^ for the change in Stokes parameters across a QT region 



for initially circular polarization were presented by Melrose, Robinson & Feletto (1995) 



Analytically, one may integrate the matrix equation Q to construct a Mueller matrix, 
which formally solves the transfer problem (Melrose 1992 Melrose fc Judge]|2004 ). Here the 
specific interest is in the position angle ?p, and it is more convenient to use ^ to derive a 
differential equation for if). 
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3.2. Natural wave modes 

The properties of natural wave modes are implicit in the square matrix in Q, and 
can be derived in terms of the eigenvalues and eigenfunctions of this matrix. The two 
non-trivial eigenvalues of the square matrix in (jij are ±zp, p = (pg + pj; + Py) 1 ^ 2 and the 
eigenvectors correspond to the two natural modes. The polarization of the natural modes 
may be described by the axial ratio, T = ay\T\, of the o mode, and the orientation ip B of 



its major axis. For the x mode the axial ratio is — 1/T at ip B + 7T. One has (Melrose & 



McPhedran 1991) 



Pq = —Ak cos2xs cos2 , 0b, Pu — —Ak cos2xs sin2-0B, pv = —Ak sm2x B , (5) 
with 

T 2 — 1 2T 
Ak = (pl + pl + p^) 1 / 2 , cos2xb = ^—^, s i n2xB = ^^. (6) 

Approximations to Ak in terms of the plasma parameters are 

ttvIvb J |cos0|, |cos^| > v B /2v, 
Ak = { (7) 

v c { u B /2v, \cos9\<u B /2u. 
The case | cos6\ < v B j2v applies within a QT region, where Ak is independent of 9. 



3.3. A conserved quantity 

As already noted, pqQ + puXJ + pyV is conserved by (|4]). In the application to Faraday 
rotation, the initial condition corresponds to linearly polarized radiation, implying V = 0, 
and circularly polarized modes, implying Pq = pu = 0, and hence PqQ + pjjU + pyV = 0. 
It is convenient to write (Q,U, V) = /(cos 2% cos 2tp, cos 2\ sin 2ip, sin 2y)- The conserved 
quantity for Faraday rotation becomes 

cos 2\ B cos 2\ cos 2(ip — ip B ) + sin 2\ B sin 2\ = 0. (8) 
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Within a QT region, the modes are elliptically polarized, and the radiation becomes 
elliptically polarized, implying that both terms in (J8|) are nonzero, so that ^ constrains 
the angle, ip — ips, between the orientation of the polarization ellipse of the radiation and 
that of the o mode. At the center of the QT region, where the modes are linearly polarized, 
sin2xB = requires cos2(-^ — ip B ) = 0, implying ip — ip B = ±tt/4. 

The conserved quantity ^ has the surprising implication that one is not free to 
specify the position angle, ip- m , of radiation incident on a QT region. One may interpret 
this as follows. The angle ipg is indeterminate in the limit of circular polarization, and as a 
QT region is entered and the modes become significantly elliptical, the angle ip B becomes 
relevant in determining the orientation of the polarization ellipses of the two natural modes. 
The requirement that the radiation be an equal mixture of the two modes determines the 
angle ip in terms of the ip B (and Xb)- Hence, ip is predetermined by the assumption that 
the radiation is an equal mixture of the two modes, and there is no freedom to specify 



3.4. Faraday angle within a QT region 

One may use Mj) to derive coupled equations that describe the changing ellipticity, 
described by Xi an d position angle, ip, within a QT region. The ellipticity changes according 
to 



dx _ i 
d£ 



2 pcos2xB sin2(?/> - ip B ). (9) 



The position angle changes according to 



cos 2 2x 



2 p[sin 2xb - cos 2xb tan 2x cos 2(ip - if) B )] = (10) 



where the second form follows by using pi). 
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Fig. 1. — An arbitrary polarization, represented in terms of Stokes parameters by q = Q/I, 
u = U/I, v = V/I is represent on the Poincare sphere by colatitute 1\ and longitude 2ip, 



with if) describing the plane of linear polarization. [After Melrose & McPhedran (1991) 



3.5. Model for a QT region 



A model for a QT region needs to specify how the plasma parameters vary with 
distance, £, along the ray path. Relevant parameters are the angles that describe the 
orientation of the magnetic field with respect to the ray path. One has cos# = at the 
center of the QT region, at £ = l\ say, corresponding to svo.2xb = 0, and the model needs to 
specify how cos# and hence sin2x# vary on either side of this zero. In a symmetric model 
sin2x B is an odd function of £ — £\. A simple example of a model is xb = Xb(£ ~ £\) with 
Xb constant near £ — 1\. An untwisted magnetic field corresponds to ipB = constant. For 
a twisted magnetic field the model must also specify how tpB varies as a function of £. An 
analogous simple model is ipB = ^b{£ ~ £\) with ip B constant near £ — l\. Mode coupling is 
strong when either xb or ipB is large compared with A A;, and mode coupling is weak when 
both are small compared with Ak. 
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3.6. Motion on the Poincare sphere 

Interpretation of the equations ^ and ( 10 ) that determine the evolution of x an d ip 
is facilitated by a pictorial description of the evolution of the polarization on the Poincare 
sphere, as illustrated in Fig. [TJ An arbitrary polarization is represented by a point on 
the sphere, with the two circular polarizations corresponding to the two poles, and linear 
polarization to the equator. The angle 2ip corresponds to longitude and the angle 2\ 
to colatitude. The orthogonal polarizations of the two modes correspond to antipodal 
points, and define a "mode axis" through the sphere, with end points at 2xb,2iPb and 
7r — 2xb,2^>b + The evolution described by Q or Q and (10) corresponds to the 



polarization point rotating relative to the mode axis. Outside QT regions, the modes are 
circular, and the mode axis is either north-south or south-north. The direction of the mode 
axis reverses each time a QT region is crossed. In the limit of strong mode coupling the 
polarization point is approximately stationary during the reversal of the mode axis. In the 
SA limit the reversal of the mode axis is slow, with the polarization point rotating rapidly 
about it many times as it moves from, say, north-south to south-north. 



The interpretation of the conserved quantity (|8j) is that the polarization point rotates 
in a great circle about the mode axis. This constraint applies specifically only to Faraday 
rotation in the SA limit, and expresses the requirement that the radiation remain an equal 
mixture of the two modes. 

In discussing Faraday rotation, the quantity of interest is ?/>, which varies periodically 
at a rate Ak. In a QL region, the great circle is around the equator, and the polarization 
point moves at constant velocity, Ak, around it. As the QT region is entered, the mode axis 



moves away from south-north, and the great circle tilts. According to (10), the variation 
in ip continues with the same period, but with a non-uniform velocity, being slower near 
the points where the great circle crosses the equator and faster where the great circle is 
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furthest from the equator. This non- uniform dif)/d£ reflects uniform circular motion around 
the great circle projected onto the equatorial plane, where the projection is an ellipse with 
angle 2ip. As the mode axis approaches the equator, the projection of the great circle onto 
the equatorial plane becomes highly eccentric. 

At the center of the QT region the mode axis is in the equatorial plane, and the 



great circle passes through both poles. According to (10), the rate of change of ip becomes 
formally infinite for sin 2% = 0, corresponding to the poles. This has a simple interpretation. 
The motion around the great circle is at constant longitude, which changes abruptly by ±7r 
as a pole is crossed: specifically, 2ip is at one value as the point moves from the south pole 
to the north pole in the western hemisphere, and changes abruptly by — 7r as it crosses the 
north pole and enters the eastern hemisphere, where it remains constant until the point 
reaches the south pole. With the projection of the great circle reducing to a line, the point 
is stationary at one end point for half the period, moving at infinite speed to the other end 
point as a pole is crossed. 

On crossing the center of the QT region to the other side, the sense of the motion 
around the ellipse formed by projecting the great circle onto the equatorial plane reverses. 
If the QT region is symmetric about its center, the increase in if> on one side of the center 
is balanced exactly by the decrease in ip on the other side. Thus, the net effect on tp of a 
symmetric (untwisted) QT region is nil. The formula (TXT) for Faraday rotation then applies 
without modification in the SA limit. A net change Aip is possible for a QT region only if it 
is asymmetric in some sense. For an untwisted field, this requires that cos 9 not be strictly 
an odd function of £ — 1\. The actual value of Aip then depends on the specific form of this 
asymmetry. 

If the magnetic field is twisted, there is a contribution to Aip, in the SA limit, due to 
the twist. This is due to the radiation remaining in the natural modes as the orientation of 
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their polarization ellipses follows the changes in longitude, 2^s, + 7r, on the Poincare 
sphere. The contribution to Aip can be estimated as ^gAI, where At is the length of 
the ray path corresponding to the QT region. Although a detailed calculation based on a 
specific model is needed to determine A£, an approximate estimate is the distance between 
the points where tan 2ip B is equal to ±1. 

In summary, the periodic variation of ip continues through the QT region, and reverses 
sense at the center of the QT region, where the requirement that there be an equal mixture 
of two modes constrains the great circle to pass through both poles, with ip having one 
of two possible values, corresponding to the longitude of the great circle in the eastern 
and western hemispheres. In a symmetric QT region, the net effect on %p is zero, implying 
that no change in ([I]) occurs in passing from the limit of strong to the limit of weak mode 
coupling. Small changes, Ai/j, can result from an asymmetry in the QT region, including a 
twist in the magnetic field. 

3.7. Polarization "reflection" at a QT region 

The foregoing discussion implies that BB's suggested "reflection" tp — > —tp at a QT 
region in the SA limit does not occur. The model used by BB (their Appendix B) assumes 
radiation incident on the QT region in an equal mixture of two modes, which I write as 
e lA( t>/ 2 e + e~* A</i//2 e x , with Acf) the phase difference between the two modes. Suppose o = R, 
x = L on the incident side, and that the components remain in their initial modes. The 
handedness of each mode reverses across the QT region, implying o = L, x = R for the 
emerging radiation. With 2ip = A<fi on the incident side, BB interpreted the interchange of 
the coefficients of the R and L polarization vectors as implying the putative polarization 
reflection, ip — > —ip. However, my analysis implies that the assumption of an equal mixture 
of the two modes requires ip — ips — ±7r/4 at the center of the QT region, effectively 
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constraining it (at the center of the QT region) to what BB call "the x = y-axis" about 
which their putative reflection occurs. There is no "reflection" : ip changes systematically in 
one sense until the center of the QT region is crossed, and then it changes systematically in 
the opposite sense. 



4. Circular polarization induced by a QT region 

In this section I point out that a QT region can lead to a partial conversion of linear 
into circular polarization. The inverse is well known: circularly polarized radiation incident 
on a QT region in the intermediate range where mode coupling is neither strong nor weak 



leads to emerging radiation that is partially linearly polarized ( Zheleznyakov & Zlotnik 
1961. 



4.1. Mode coupling across a symmetric QT region 

Mode coupling is determined by a coupling coefficient, Q x , that is effectively the ratio 
of two rates: that at which the shape or orientation of the polarization ellipse of the natural 
modes changes, and the rate at which the modes get out of phase. Mode coupling is strong 
when this ratio is large and weak when it is small. In QL regions, where the polarization of 
the natural modes does not change, mode coupling is weak. As a QT region is approached, 
the polarization becomes significantly elliptical and Ak decreases: Q x increases due to both 
effects. In the case of weak coupling one has Q x <C 1 throughout a QT region. Consider 
intermediate mode coupling, with Q x ^> 1 at the center of the QT region, and Q x > 1 for 
some distance on either side of the center. In this case, weak mode coupling applies to 
radiation approaching the QT, with Q x increasing from <C 1, and passing through unity, say 
at £ = £-. Mode coupling is then strong across a region, £-<£<£+ say, that includes the 
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center, £ = £i, of the QT region, becoming weak when Q x passes through unity at £ = £ + . 

Consider an idealized model in which the components in the two modes follow the 
changing polarization of the modes up to £ = £-, at which point the polarization of the 
radiation is frozen-in, remaining unchanged across the center of the QT region and ceasing 
to be frozen-in at £ = £ + , after which the components in the two modes again follow the 
changing polarization of the modes. In this model, the condition ^ is satisfied for £ < 
Let the polarization at £ = £_ be described by \ = X- an d ip = ij)— The radiation is an 
equal mixture of the two modes so that ^ is satisfied with \ — X—> V* = ' t P- an d Xb = Xb-, 
ipB = ipB~- The model implies that x = X- an d ip = ip- remain constant for < £ < £ + . 

In terms of the vector notation in which Q becomes dS/d£ = p x S, the invariant 
p ■ S/p applies for £ < £_ and for £ > £ + . The constant quantity over the range £_<£<£+ 
is S = S_, with p varying and p = p± different at £ — £±. 

For a symmetric QT region, at £ = £ + one has xb+ = ~Xb- an d ipB+ = V's- + tt, 
corresponding to p + = — p . The condition ^ is then satisfied at £ = £ + . In the vector 
notation, one has p + ■ S = —p ■ S_ = at £ = £ + . In a symmetric model, the decrease in ip 
for £ > £i is a mirror image of the increase for £ < l\. The changes across a symmetric QT 
region cancel exactly. This implies that the standard formula ([T]) for RM applies without 
modification irrespective of the strength of the mode coupling, provided that the QT region 
is symmetric. 

4.2. Effect of an asymmetry 



Suppose that the QT region is asymmetric in the sense that at £ = £ + one has 
Xb+ 7^ ~Xb- and/or ip B+ ^ ip B _ + n, that is, p + ^ —p in the vector notation. This is 
obviously the case if there is a twist in the magnetic field, when ipB+ differs from i^b- + tt 
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by the angle through which the magnetic field twists between £ = £_ and £ = £ + . There can 
also be an asymmetry causing xb+ to differ from — \b- due to cos^ not being a strictly 
odd function of £ — L . 



In the absence of symmetry, the condition (|8j) is not satisfied at £ = £ + . To be specific, 
suppose one has 

cos 2\b+ cos 2x- cos 2(ip — ipB+) + sin 2\b+ sin 2%_ = 5 + (11) 
at £ = £ + , where \ = X-i ^ = i ) - applies to the frozen- in polarization which satisfies ((8| 



with xb = Xb~, 4>b = 4>b~- As ips, Xb change with £ > £ + , \ an d ip vary such that (11) 
continues to be satisfied with fixed 5 + . In terms of the vector notation, for £ > £ + one has 
p ■ S/p constant at its value at £ = £ + . 



The implication of (11) is that the radiation is an unequal mixture of the two modes 
for £ > £ + . For £ ^> £ + , in the subsequent QL region, this implies that the radiation is 
an unequal mixture of two opposite circular polarizations, implying a degree of circular 
polarization 5 + . The handedness of the resulting polarization depends on the sign cry and 
the nature of the asymmetry, e.g., the sense of twist of the magnetic field. One expects it 
to be random: for a large number of QT regions the signs 5 + > and 5 + < would occur 
with equal probability. 

A detailed model is required to estimate how the magnitude of 5 + varies with frequency. 
The following arguments suggest that 5 + is a maximum near the transition frequency, 
vt- Mode coupling becomes stronger with increasing v > ut, implying that for v ^> 
the freezing of the polarization begins in the QT region when the modes are only slightly 
non-circular. This corresponds to small cosxb-, sin%_ and hence small S + , with 5 + 
decreasing as the mode coupling becomes stronger. At low frequencies, v vt, there 
is negligible mode coupling, and negligible freezing-in implies negligible induced circular 
polarization. The optimum conditions for induced circular polarization occur for v w z/^, 
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where cosxb- and sinx_ are not small, and the distance, Ai = £ + — is large enough for 
a twist in the magnetic field to have a significant effect on the orientation of the polarization 
ellipses of the natural modes. 

5. Observable consequences 

I comment on the possibility of observing three effects: the change Aip, circular 
polarization due to mode coupling, and the quantity defined by §2§. 

5.1. Effect of A^p y£ 

For a change, due to mode coupling changing from strong to weak, to be observable, 
the transition frequency, vt, must be in the range of observation. BB estimated Vt for the 
ISM, and argued that under favorable circumstances it may be in the range of observational 
interest, 100MHz-l GHz. Let me assume that this is the case, and consider what would be 
observed. 

If BB's suggestion were correct that RM is given by ^ rather than in the SA 
regime, the change in RM with decreasing frequency would be potentially very large, with 
([TJ applying at high frequency and ^ at low frequency. However, ^ is not relevant 
to Faraday rotation, and any actual change in RM from high to low frequency is much 
smaller than (|2| would suggest. A change Aip due to a twist in the field can be written 
as Aip = ipsAi, where A£ oc A is the thickness of the QT region. This effect can be 
observed by determining RM at short wavelengths, and using it to identify the position 
angle ip at the source (by taking the limit A — > 0), and then repeating the procedure at 
longer and longer wavelengths. The prediction is that one would find that the inferred 
ipo has a contribution oc ipB^ that depends on the wavelength band used to determine it. 
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Observation of such a change would enable one to determine the transition frequency, Ut, 




A complication is that one expects there to be many QT regions along the line of sight. 
Statistically, half of the changes would have Aip > and half have Aip < 0. Hence, if the 
number of QT regions is large, one expects their cumulative effect to sum to zero, leaving 
the standard result essentially unchanged for v <C ut- 



5.2. Circular polarization due to mode coupling 

Circular polarization induced by a QT region on incident linearly polarized radiation 
is a new effect that has not been recognized previously It has some similarity to 



scintillation-induced circular polarization (Macquart & Melrose 2000), but is not the same 



effect. There is a partial analogy with linear polarization induced by a QT region on 



incident circular polarization (Zheleznyakov & Zlotnik 1964). A notable feature of the 



analogy is that both effects are maximized when mode coupling is intermediate between the 
strong and weak limits. This maximum should be near the transition frequency, v ~ v?. 
Observationally, the effect should be confined to a bandwidth of order vt about Vt for each 
QT region. 

For a QT region to induce circular polarization it must be asymmetric, and a plausible 
cause for the asymmetry is a twist in the magnetic field. One would expect this effect to be 
observable when there is only one QT region along the ray path. In practice one expects 
many QT regions, and for a statistically large number of QT regions, the degree of circular 
polarization from each is random, in the sense that there is no systematic sign, and hence 
would tend to average to zero. 
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5.3. Comparison with solar application 

There is a long-standing problem relating to bipolar regions in the solar corona, where 
x mode radiation from opposite footpoints must be oppositely circularly polarized on 
emission. The radiation from the more distant footpoint must pass through an extra QT 
region compared to radiation from the nearer footpoint. A prediction is that for v < vt, the 
two footpoints should have the same handedness. Estimates of vt for propagation through 
the solar corona suggest sufficiently high values that one should observe this change above 



about 100 MHz. However, the effect is not observed (White, Thejappa & Kundu 1992), 
suggesting that the value of is overestimated. 

In my opinion, the likely resolution of this long-standing problem is that the magnetic 
field has a small-scale random component, and that what is identified as a single QT region 
where the mean field passes through 9 = n/2 consists of a large number of mini-QT regions 
where the local value of 9 passes through tt/2. The strength of mode coupling depends on 
L in (|3j) , and this is very much smaller for the mini-QT regions than would be estimated 
for the mean field. A random small-scale component of the magnetic field in the ISM would 
have an analogous effect, reducing the likely value of to much less than the estimated 
range 100 MHz- 1 GHz of observational interest. 

Suppose the magnetic field has a systematic twist, as well as a random component. 
The additional rotation, Aip, and the induced circular polarization, 5 + , each have opposite 
signs at consecutive mini-QT regions. The number of mini-QT regions is necessarily odd, 
and the net effect should be of the same order as for a single mini-QT region. This suggests 
that inclusion of a random component makes it more difficult to observe Aifj, which requires 
v <^vti but has less effect on observing 5 + , which requires v ut- 
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5.4. Polarization-dependent DM 

Formula ^ is not relevant to Faraday rotation, but the quantity ^ can be measured 
in principle through a relative time delay between the o mode and x mode components. 
If the SA limit applies at every QT region between the source and the observer, the 
components in the two modes propagate independently, and (|2| determines the accumulated 
phase difference between them. The difference in group velocity between the o and x mode 
components is proportional to n D — n x , and this leads to a time delay in the arrival of the 
oppositely circularly polarized component of a pulse. This relative delay between oppositely 
polarized components increases oc A 3 , with the coefficient cx RMsa- Observationally, 
this implies a correction of opposite sign to the DM for oppositely circularly polarized 
components in a pulse. Although observing a relative time delay between oppositely 
circularly polarized components is possible in principle, it is essential to have a well-defined 
pulse to detect. A potential complication is dispersive smearing, which tends to smear out 
a pulse, with this smearing increasing as the frequency decreases. 



6. Conclusions 

This investigation of the effect of mode coupling at a QT region on Faraday rotation 
was motivated by a recent discussion by BB, who argued that at low frequencies, where 
mode coupling is weak, the standard formula (IT]) for Faraday rotation is replaced by ([2]). 
This is based on a misinterpretation of what happens at a QT region; ^ is not relevant to 
Faraday rotation. The analysis given here implies that the standard formula applies 
irrespective of the strength of mode coupling at QT regions, provided that each QT region 
is symmetric. An asymmetry can induce a small correction, Aip, that needs to be added 
to RMA 2 in the SA limit. The sign of this correction is essentially random, and the net 
correction over many QT regions should average to zero. 
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There are three other new results to emerge from the present investigation, one that is 
theoretical, and two of possible observational relevance. 

The theoretical result relates to the limit of weak mode coupling (BB's SA limit), when 



new result is that this invariant is equal to zero for Faraday rotation; this particular value 
corresponds to an equal mixture of the two modes. One implication is that the polarization 
point on the Poincare sphere moves around a great circle, whose orientation changes as the 
polarization of the natural modes changes. A particular implication follows by considering 
this motion at the center of the QT region where the natural modes are linear: the great 
circle then passes through both poles, corresponding to cyclic conversion of linear to circular 
polarization as in a quarter- wave plate. The implication is that the position angle, if), at 
the center of the QT region is predetermined by this invariant and the orientation of the 
magnetic field at that point. This invalidates the simple model used by BB to justify ^ 
for Faraday rotation because this model pre-supposes that the initial orientation of the 
plane of polarization is arbitrary. More specifically, the reflection, ip — > —if>, about the 
u x = y-axis" that BB argued occurs at a QT region is a misinterpretation: the u x = y-axis" 
is the position angle of the great circle at the center of the QT region. 

The development of circular polarization at a QT region can be understood in terms of 
this invariant. At intermediate frequencies, v « vt, mode coupling is weak on the edges of a 
QT region, and strong in a region about its center. The quantity (pqQ + PuU + pyV)/AkI 
is invariant only when mode coupling is weak. When mode coupling is strong, across a 
region about the center of the QT region, the polarization of the radiation is frozen-in. Any 
asymmetry across this region, such as a twist in the magnetic field, implies that the values 
of (pqQ + puU + pyV)/Ak I at the two points where the transition occurs between weak 
and strong mode coupling are not equal. An initial value of zero implies a non-zero final 



the polarization transfer equation Ml) implies an invariant 
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value, and hence an unequal final mixture of the two modes. On emerging from the QT 
region this becomes an unequal mixture of two circular components, giving a net degree of 
circular polarization, S+. This induced circular polarization should be a maximum near the 
transition frequency, vt- 

Although observation of the effects described here by Aip and 5 + encounter difficulties, 
these may be less serious for 5+, which requires v « i/?, than for Aip, which requires v <C Ut- 
The signs of both and 5 + depend on details that suggest that they should average to 
zero over a larger number of QT regions. Realistically, the effects could be observed only if 
at most a few QT regions contribute. 

The quantity defined by the alternative formula (|2| proposed by BB for RM, although 
not relevant to Faraday rotation, is observable in principle as a relative time delay between 
the o mode and x mode components in the SA limit. Observationally, this effect implies 
a correction of opposite sign to the DM for oppositely circularly polarized components in 
a pulse. Any measurement of this effect would give intrinsically new information of the 
properties of the ISM. Because this effect favors low frequencies, dispersive smearing is 
likely to be a problem in any attempt to measure it for pulsars. 

I thank Alex Judge for helpful discussions. 
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